A DISCRETE GAUSS-GREEN IDENTITY FOR UNBOUNDED LAPLACE 
OPERATORS, AND THE TRANSIENCE OF RANDOM WALKS 



PALLE E. T. JORGENSEN AND ERIN P. J. PEARSE 



Abstract. A resistance network is a connected graph (G, c). The conductance function c xy 
weights the edges, which are then interpreted as resistors of possibly varying strengths. The 
relationship between the natural Dirichlet form £ and the discrete Laplace operator A on a 
finite network is given by £(u, v) = (u, Av)2, where the latter is the usual C 2 inner product. 
We extend this formula to infinite networks, where a new (boundary) term appears. The 
Laplace operator is typically unbounded in this context; we construct a reproducing kernel 
for the space of functions of finite energy which allows us to specify a dense domain for 
A and give several criteria for the transience of the random walk on the network. The 
extended Gauss-Green identity and the reproducing kernel also yield a boundary integral 
representation for harmonic functions of finite energy; this representation is the focus of a 
forthcoming paper. 
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1. Introduction 

Our purpose is to extend to infinite weighted graphs the Gauss-Green identity 

I VmVvc/V = - I uAvdV+ | uf-dS, (1.1) 
Ja Jn Jen 11 

which expresses a fundamental relationship between analytical and geometrical notions 
of boundary. When the domain Q is a network (connected weighted graph) G, then the 
boundary bdG = dQ. will be a limit object similar to the Martin boundary of a Markov 
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chain (in particular, the boundary is not some prescribed subset of the domain). In this 
discrete context, the left side of (1.1) corresponds to the (Dirichlet) energy form 

£(«, v) :=- ^ J] ^("W - «(y))(v(*> - v(y)), (L2) 

.tec veG° 

where G° is the vertex set of G. Furthermore, the Laplacian on the right side of (1.1) is 
given pointwise by (Av)(*) := Yiy~x Cxy(v(x) - v(y)), where x ~ y means that there is an 
edge between x and y, and c„ = c yx > is the weight ("conductance") of this edge. For 
finite networks (or finite subnetworks of infinite networks), an elementary computation 
shows that the energy S and Laplacian A are related by the formula 

S(u, v) = Y "Av, (1.3) 

G° 

so that the boundary term corresponding to J gn uj^ dS in (1.1) vanishes. However, this 
is not always the case for infinite networks. Let us define a function on G° to be har- 
monic iff Au(x) = for every x e G°. Then it is clear that formula (1.3) cannot hold 
for networks which support nonconstant harmonic functions of finite energy (also called 
harmonic Dirichlet functions), because this would imply 

< &{h, h) = Y h(x)Ah(x) = Y /! W -0 = 0, (1.4) 

xeG° xeC° 

where the strict inequality for nonconstant functions follows from the connectedness of the 
network and the fact that (1.2) is a sum of nonnegative terms for u = v = h. 

The main result of this paper is Theorem 3.12, in which we resolve (1.4) and describe 
precise conditions for an infinite network to have a nonvanishing boundary term. In partic- 
ular, with domfi = {u • &(u) < oo}, we have 

S(u,v) = Y u ^ v + Y U ^' for all u e domfi, v e M, (1.5) 

G° bdG 

where M is a certain domain for the Laplacian, and bd G and are explained in Defini- 
tion 3.9 and Definition 3.10. In particular, it should be noted that bdG is not a subset of 
G° or some larger ambient graph; instead, this object may be understood as a collection 
of certain equivalence classes of geodesies in the spirit of [Car72, Car73]. It is the focus 
of [JP09a] to understand bdG as a measure space; see also [JP09d]. 

Until now, results in the literature on infinite networks have focused on finding suffi- 
cient conditions to ensure the boundary term vanishes; see the discussion of the literature 
just below. This paper also provides several conditions under which this occurs (see Re- 
mark 3.16), but more importantly, it also provides conditions under which the boundary 
term is convergent and nonzero. The correct determination of M is actually the primary 
challenge for this; see Definition 3.4. We additionally give methods to compute the value 
of such a boundary term. 

Theorem 3.12 implies a boundary representation for harmonic functions of finite en- 
ergy: Corollary 3.13 states that for a harmonic function u with &(u) < oo, one can recover 
the values of u via 
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m « = Z m It +m(o) - (L6) 

bdG 

For now, we consider bdG and ^- only as objects defined in terms of limits. However, 
Corollary 3. 1 3 is the first step in the development of a boundary theory for harmonic func- 
tions of finite energy (akin to the Martin boundary theory for nonnegative harmonic func- 
tions, or the Poisson boundary theory for bounded harmonic functions). The paper [JP09a] 
formalizes the construction of bdG as a measure space (see 3.11). The resulting object 
bdG serves as a support set for a measure which represents the harmonic functions of 
finite energy, in the sense that 

u(x) = I k(x,d/3), 

JbdG 

for any harmonic function u of finite energy, where k(x, dp) is a kernel whose second 
argument is a measure supported on bdG. Consequently, we expect that the resistance 
boundary bd G will be useful for studying induced Dirichlet forms on the boundary of an 
infinite resistance network. If tp is a function on bd G for which the harmonic extension 
Htp to G has finite energy, then 



SbdG(v) := £>(H<p) 

defines the induced Dirichlet form on the resistance boundary of G; see [Kig09a]. This 
provides an approach to doing analysis on spaces which are otherwise difficult to deal with, 
including self-similar fractals, Julia sets, and limit sets of Kleinian or Fuchsian groups. 
See [DS01, DS99, Kai03, Kig09a, LW09, JLW, RT10]. 

Remark 1.1 (Motivational examples). The boundary term does not vanish for homoge- 
neous trees [Car72, Car73, Kig09a]; tesselations of the hyperbolic disk and graphs which 
are roughly isometric to hyperbolic spaces [LP09, §2 and §9]; graphs which are uniformly 
embedded in the hyprbolic disk [CW92]; and contacts graphs of circle packings in the 
plane [BS96]. Some examples are discussed in §6. 

As an elementary example where the boundary term does not vanish, consider the fol- 
lowing network discussed in Example 6.3: let G be the network whose vertices are the 
integers, where n is connected to « - 1 by an edge of conductance c n>n -\ - c" for some 
c > and for n > 0, and symmetrically for « < 0: 



Lemma 6.5 gives a nonconstant harmonic function of finite energy h, and Remark 6.6 
shows that for this function, the boundary term is 2 bdG hjj^ = 1. 

For c = 2, this example is closely related to the example of the binary tree network 
(with Cxy = 1 for all x,y e G°). However, the network (1.7) provides an example where 
the Laplacian is unbounded and in fact is not self-adjoint. While the Laplacian as defined 
above is always a Hermitian operator, it is not even essentially self-adjoint in this example; 
that is, there exist multiple nontrivial self-adjoint extensions. 

In [JP09e, §] and [JP09c, Ex. 13.39], we show that the network (1.7) supports a defect 
vector, that is, a function / satisfying S(f) < oo and Af(x) = -f(x). (It turns out this 
function is also bounded.) For this function, 
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J] f(x) A f(x) = £ fix)i-fix)) = - J] l/MI 2 (1.8) 

A-eG° ,i-eC° -veC 

and one would have a counterexample. This highlights the importance of M; note that 
/ £ M. Consequently, a key point of Theorem 3.12 is that for v € M, the two sums in 
(1.5) are both finite. 

Notes on the literature. Random walks on graphs (and trees in particular) comprise an 
old and well-studied subject and we will not attempt to give complete references. However, 
we recommend [DS84,Doy88,LP09,LPW08] for introductory material and [Lyo83,Car73, 
WoeOO], and the foundational paper [NW59] for more specific background. With regard to 
infinite graphs and finite-energy functions, see [Soa94,SW91,CW92,Dod06,PW90,PW88, 
Woe86, Tho90, BS96]. Effective resistance and resistance metric are studied extensively 
in [Kig03,Kig09b], and also in [KigOl]. 

For finite (sub) networks, a formula similar to (1.5) appears in [DK88, Prop 1.3], and 
another similar formula appears in [RSV97]. However, these authors do not give an exten- 
sion of this formula to infinite networks. 

Infinite networks are discussed in [KY89, Thm. 4.1], where the authors give some con- 
ditions under which (1.3) extends to infinite networks. However, these authors do not 
discuss the case when (1.3) must be replaced by (1.5). The emphasis of the present paper 
is on this latter case because the (nonvanishing) boundary term leads to a boundary repre- 
sentation for harmonic functions of finite energy, as well as detailed information about the 
structure of the space of functions of finite energy, and several computational techniques 
discussed below. Furthermore, the conditions provided by Kayano and Yamasaki may be 
highly nontrivial to verify. By contrast, we develop a dense subspace of functions to which 
the formula may be applies, and in [JP09b], we show that these functions are relatively 
easy to compute. 

A reproducing kernel similar to energy kernel {v x } xeC o constructed in §2.1 appears 
in [MYY94] under the name Kuramochi kernel. Indeed, the Kuramochi kernel can be 
recovered as a special case of the present work. As an advantage of the present approach, 
we note that our formulation puts the Green function in the same space as the reproducing 
kernel elements v x . Also, the geometry of the Hilbert space studied in this paper is much 
simpler, in comparison to that induced by the inner product studied in [MYY94]; see §4. 1 . 

It is easy to think of the formula (1.3) as <5(m, v) = (u, Av)p, and hence the reader may 
be tempted to see our results as part of Kato's theory. Recall that [Kat95] gives a bijective 
correspondence between the family of closed symmetric quadratic forms fiona Hilbert 
space 'H and the family of nonnegative definite self-adjoint operators A on *H via 

dom£ = dom( VZ), 8 (m, v) = ( VAm, VAv)^ . 

See [FOT94, Thm. 1.3.1] and its corollary. However, if one requires u,v e ( 2 iG°), then 
one does not usually obtain a dense subspace of the functions of finite energy. On the 
other hand, requiring that u and v have finite energy does not guarantee that either is in 
£ 2 iG°). Furthermore, even if Kato's theory could be adapted to accommodate this issue of 
domains, it still does not account for the second term on the right side of (1.5). 

The reader may also notice a similarity between the Hilbert space 'Hg which is the focus 
of the present study, and the extended Dirichlet space T e studied in [FOT94, § 1 .5— § 1 .6] 
and elsewhere. However, the resemblance is only superficial and in general there is no 
relation between the two. Note that T e is a Hilbert space if and only if the Dirichlet form 
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is transient, and that the ( 2 functions are dense in T e [FOT94, Thm. 1.5.3]. These are both 
false for 

Applications to mathematical physics. Due to the discrete nature of (1.5), this formula 
has applications beyond the problems of geometry. In particular, there is an immediate 
interpretation in terms of functions on X x X (where X is a discrete set) which is naturally 
understood in terms of matrices; see [JP09e]. 
The symmetric form of the identity (1.1) is 

f(uAv-vAu)dV = f (u§±-v^)dS, 
Jo. Jan 

and in a Hilbert space, this takes the form 

(<A*w,v>-<w,A*v>) =j8(z*,v), w.vedomA*, (1.9) 

where A* is the adjoint of A with respect to the inner product of f H, and the right-hand 
side of (1.9) defines a boundary form in the sense of [DS88]. Many problems for trans- 
formations in linear algebra may be formulated with the use of symmetric matrices and 
the determination of [5 in (1.9) is an immediate question. In finite dimensions, all inner 
products are equivalent, and the boundary term is trivial. However, this is not so for prob- 
lems from statistics, particle models, or infinite weighted graphs which require analysis in 
infinite dimensions [ADV09,Han96]. These analyses entail a choice of inner product, and 
therefore by completion, a Hilbert space Iri, and rigorous study of boundary problems for 
an infinite discrete system X typically dictates a particular Hilbert space. If X is the vertex 
set of an infinite graph, it is tempting to take H = ( 2 {X) as the preferred Hilbert space, but 
it turns out that £ 2 {X) gives a trivial boundary theory for X, in the sense of (1.5). 

In this context, results of the present paper may be applied to certain discrete analogues 
of Laplace (or Schrodinger) operators. However, our framework applies to a variety of 
other discrete models as well, for example to operators in statistical physics generating the 
dynamics of particle spin-flips; see e.g., [CL07,FHM03,Lig78]. For example, consider the 
following close relative of the Laplace operator which plays a key role in spin dynamics: 

(W)(J7) = Yj WW) - /to)). n 6 {°> 1 ) A 

where r\ is a configuration of particle spins, 

w s f i - y = x, 

\T](y), yl=x, 

and c{x, rf) is a nonnegative function translation-invariant which is invariant on the lattice 
A. The case A = Z is discussed in [Lig78]. 

Outline. The rest of this paper is as follows: 

§2 develops the Hilbert space structure of the functions of finite energy, in terms of a 
certain reproducing kernel {v A } for dom£. This kernel is necessary for formulating the 
hypotheses of Theorem 3.12, but it also gives insight into the Hilbert space structure of 
dom£. Additionally, this kernel (and the associated framework developed in the present 
paper) is used to study effective resistance metric on infinite networks in the forthcoming 
paper [JP09b], to construct a new boundary integral representation for harmonic functions 
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of finite energy in the forthcoming paper [JP09a], and to obtain information about the 
spectrum of the Laplacian in this context. 

§3 contains the complete statement and proof of Theorem 3.12, as well as the precise 
definitions required. In particular, we discuss the class of monopoles, which can be thought 
of as an extension of Green's kernel, and their relationship to the reproducing kernel(s) 
introduced in the previous section. Corollary 3.13 introduces a boundary representation 
for harmonic functions studied further in [JP09a]. 

§4 relates the boundary term £bdc u jti to tne transience of the network. We indicate 
how the characterizations of transience due to [Lyo83] and [NW59] appear in the present 
context, and introduce three new equivalent characterizations of transience, stated in terms 
of our main formula Theorem 3.12 and the operator-theoretic properties of A. 

§5 concerns implications and applications of (1.5), and some facts about useful special 
cases. We recover several known facts about finitely-supported functions and harmonic 
functions; the use of the energy kernel {v A } and other techniques introduced in previous 
sections allow for shorter and easier proofs. We also give a result indicating how one can 
extend Theorem 3.12 to networks which contain vertices of infinite degree. 

§6 gives a family of examples that illustrates the properties we discuss, and several 
of our results. We prove that these simple networks support monopoles and harmonic 
functions, and discuss why A may not be self-adjoint on these networks for appropriate 
choices of c. 

1.1. Basic terms. We now proceed to introduce the key notions used throughout this pa- 
per: resistance networks, the energy form <5, the Laplace operator A, and the elementary 
relations amongst them. 

Definition 1.2. An resistance network is a connected graph (G, c), where G is a graph 
with vertex set G°, and c is the conductance function which defines adjacency by x ~ y 
iff Cxy > 0, for x,y e G°. We assume c xv = c yx e [0, oo), and write c(x) := Y* y ~x c xy We 
require c(x) < oo, but c(x) need not be a bounded function on G . The notation c may be 
used to indicate the multiplication operator {cv)(x) := c{x)v(x). 

In this definition, connected means simply that for any x,y e G°, there is a finite se- 
quence with x — xo, y — x n , and c x ._ lX . > 0, i = l,...,n. Conductance is the 
reciprocal of resistance, so one can think of (G, c) as a network of nodes G° connected by 
resistors of resistance c~ v ! . We may assume there is at most one edge from x to y, as two 
conductors cjL and c XY connected in parallel can be replaced by a single conductor with 
conductance c xv = c' v + cL. Also, we assume c xx — so that no vertex has a loop, as 
electric current would never flow along a conductor connecting a node to itself. 1 

Definition 1.3. The Laplacian on G is the linear difference operator which acts on a func- 
tion v : G° -> R by 



A function v : G — > R is called harmonic iff Av = 0. The Laplacian may be written 
A = c — T, where (T v)(x) :— Yjv~x c xvv(y) is the transfer operator. 



Nonetheless, self-loops may be useful for technical considerations: one can remove the periodicity of a 
random walk by allowing self-loops. This can allow one to obtain a "lazy walk" which is ergodic, and hence 
more tractable. See, for example, [LPW08, LP09]. 




(1.10) 
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The network Laplacian defined above should not be confused with the normalized 
Laplace operator defined by c~'A, which appears frequently in the probability literature 
(e.g. [DS84]), nor the version c _1 ^ 2 Ac -1 ^ 2 which appears in the literature on spectral graph 
theory (e.g., [ChuOl]). Also note that we have adopted the physics convention (so that 
the spectrum is nonnegative) and thus our Laplacian is the negative of the one commonly 
found in the PDE literature. 

Remark 1.4. The reader may wonder why we have elected to define the Laplacian as 
above, instead of one of the other mentioned options. A primary reason for our use of 
the unnormalized Laplacian is that it may be unbounded, and therefore have a delicate 
spectral structure which is not visible in the bounded case. In particular, the Laplacian 
given by (1.10) may fail to be essentially self-adjoint. This means that A has multiple 
self-adjoint extensions; or equivalently, nontrivial defect spaces. The issue of essential 
self-adjointness is closely related to the question of whether or not the graph is stochasti- 
cally completeness, and in some contexts the two are equivalent. R. Wojciechowski has 
recently given a characterization of stochastic completeness for certain classes of graphs; 
see [Woj07, Woj09, WojlO]. Also, M. Keller and D. Lenz characterize stochastic complete- 
ness for regular Dirichlet forms on discrete sets in [KL10], and relate this to the stochastic 
completeness of graphs and subgraphs. Additionally, A. Weber has recently given a condi- 
tion for stochastic completeness of graphs which may be understood as a weak curvature 
bound [Web08]. See [JP09e] for further definitions and discussion. We expect the results 
of the present paper to be useful for studying stochastic completeness. 

Definition 1.5. An exhaustion of G is an increasing sequence of finite and connected sub- 
graphs {Gk}, so that Gk £ Gk+\ and G = [j Gk- 

Definition 1.6. The notation 

xeG° -«eG a 

is used whenever the limit is independent of the choice of exhaustion {Gk} of G. This is 
clearly justified, for example, whenever the sum has only finitely many nonzero terms, or 
is absolutely convergent as in the definition of £ just below. 

Definition 1.7. The energy of functions u, v : G° — > R is given by the (closed, bilinear) 
Dirichlet form 

£(«, v) :=\- J] Z ^("W - "WXvM - v(y)), (1.12) 

.teG° veG° 

with the energy of u given by £(m) := £(m, u). The domain of the energy is 

dom£ = {u : G° -> R-£(w) < oo). (1.13) 

Since c xy = c yx and c XY = for vertices which are not adjacent, there is exactly one 
term in the sum in (1.12) for each edge in the network when we include the factor of | to 
prevent double-counting. Note that (1.12) implies 

£(5.0 = c(jc), and 8(6 x ,6 y ) = -c^, (1.14) 

where 6 X is a (unit) Dirac mass at x e G°. Consequently, the Dirac masses S x and 6 y are 
orthogonal with respect to energy if and only if x y. 
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The following proposition may be found in [Str06, §1.3] or [KigOl, Ch. 2], for example. 

Proposition 1.8. The following properties are readily verified: 

(1) ker£ consists precisely of the constant functions. 

(2) (Polarization) &(u, v) = \[£>(u + v) - &(u — v)]. 

(3) (Markov property) S([u]) < S(u), where [u] is any contraction 2 of u. 

Proposition 1.9. Let G° be finite. Then 

S(u, v) = ^ u(x)Av(x), (1.15) 

A-eG° 

and all harmonic functions of finite energy are constant. 

Proof. The computation is elementary, but we include it for later reference: 

u( x )Av(x) = \Yj "(*) Av to + \ "Cy)Av(y) 

A-eG° xeG» yeC" 

= \ X X c xy u(x)(v(x) - v(y)) - i ^ ^ c vv «(y)(v(x) - v(y)) 

A-eG° veG •*->' 

= \ X C «(«( JC ) - u(y))(v(x) - v(y)), 

agg y-^ 

If h is harmonic, then apply this formula to get S(li) = 2 xe c° h(x)Ah(x) = 0, and the latter 
result follows from (1) of Proposition 1.8. □ 

Connectedness is implicit in the calculations behind Proposition 1.8 and Proposition 1 .9; 
recall that all resistance networks considered in this work are connected. Our main result 
is an extension of the above formula for S(u, v) to infinite networks in Theorem 3.12, 
where the formula is more complicated: (1.5) indicates the presence of a "boundary term" 
ZbdG u 7it- I* ' s snown m Theorem 4.6 that the presence of the boundary term corresponds 
to the transience of the random walk on the underlying network, that is, the Markov process 
with countable state space G° and transition probabilities p(x,y) :- c xy /c{x). 

In the traditional study of Dirichlet forms [FOT94] or more general quadratic forms 
[Kat95], one would write S{u, v) = (u, Av)p and consider the space of functions with 
||m||2 + &{u) < oo. In classical potential theory (or Sobolev theory), this would amount to 
working with the class of functions satisfying ||/'||2 < °°, but abandoning the requirement 
that II/H2 < 00. In our context, this is counterproductive: many of the most interesting 
functions in domfi are not in { 2 (G°); see Corollary 5.6. Consider that if at least two 
connected components of G\ {0} are infinite, then v x £ I 1 for vertices x in these components, 
where v x is an element of the energy kernel; see Definition 2.6. Also, in consideration of 
(1.8), it is impossible for a defect vector of A to be in t 1 . 

2. The energy Hilbert space 

In this section, we study the Hilbert space Ifs of (finite-energy) voltage functions, that 
is, equivalence classes of functions u : G a — > C where u v iff u - v is constant. On this 
space, the energy form is an inner product, and there is a natural reproducing kernel { v x } x£ qd 



"In this context, a contraction is any operation u h-» [u] such that |[«](*) - [«](y)| < \u(x) - u(y)\ whenever 
x ~ y. For example, let [h] := minj 1, max{0, u}}. 
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indexed by the vertices; see Corollary 2.7. Since we work with respect to the equivalence 
relation defined just above, most formulas are given with respect to differences of function 
values; in particular, the reproducing kernel is given in terms of differences with respect 
to some chosen "origin". Therefore, for any given resistance network, we fix a reference 
vertex o e G° to act as an origin. It will readily be seen that all results are independent 
of this choice. When working with representatives, we typically abuse notation and use u 
to denote the equivalence class of u. One natural choice is to take u so that u{o) — 0; a 
different but no less useful choice is to pick k so that v = outside a finite set as discussed 
further in Definition 2.10. 

Let 1 denote the constant function with value 1 and recall that ker£ = CI. 

Definition 2.1. The energy form & is symmetric and positive definite on dom£. Then 
dom£/Cl is a vector space with inner product and corresponding norm given by 

(u,v) s :=B(u,v) and ||w|| £ := 8(u, u) 1 ' 2 . (2.1) 
The energy Hilbert space Us is the completion of dom£/Cl with respect to (2.1). 

It can be checked directly that the above completion consists of (equivalence classes 
of) functions on G° via an isometric embedding into a larger Hilbert space as in [LP09, 
MYY94] or by a standard Fatou's lemma argument as in [Soa94]. 

Remark 2.2 (Four warnings about < H&). 

(1) 7Ye has no canonical o.n.b.; the usual candidates {S x } are not orthogonal by (1.14), 
and typically their span is not even dense; cf. Corollary 2.17. 

(2) Pointwise identities should not be confused with Hilbert space identities; see Re- 
mark 2. 13 and Lemma 3.7. 

(3) Multiplication operators are not generally Hermitian, as we show in the next 
lemma contrasts sharply with more familiar Hilbert spaces. 

(4) There is no natural interpretation of as an £ 2 -space of functions on the vertices 
G° or edges G 1 of (G, c). 

With regard to (4), it should be noted that He does contain the embedded image of { 2 {G°, /i) 
for a certain measure p, but this space is not typically dense, and almost never equal to He- 
Also, "Hg embeds isometrically into a subspace of i 2 (G l ,c), but it generally nontrivial to 
determine whether a given element of £ 2 (G l ,c) lies in this subspace. "He may also be 
understood as a I 2 space of random variables [JP09c, §15.1] or realized as a subspace of 
L 2 (S',F), where S' is a certain space of distributions [JP09a], but both of these construc- 
tions are beyond the scope of the present paper. 

Lemma 2.3. If <p : G° — > R and denotes the multiplication operator defined by 
{M,pU)(x) — tp(x)u(x), then is Hermitian if and only if — k\ for some k G R. 

Proof. Choose any representatives for u, v e 7Yg. From the formula (1.12), 

{M^u, v) 6 = r ^ c v( ( P( x M x ) v ( x ) - <p(x)u(x)v(y) - ip(y)u(y)v(x) + ip(y)u(y)v(y)). 

x,y£G° 

By comparison with the corresponding expression, this is equal to (u, M^v)g iff {<p{y) — 
tp(x))u(y)v(x) = (ip{y) - ip{x))u(x)v(y). However, since we are free to vary u and v, it must 
be the case that ip is constant and hence ip — in "Hg. The converse is trivial. □ 
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2. 1 . The evaluation operators L x and the reproducing kernel { v x ). 

Definition 2.4. For any vertex x e G°, define the linear evaluation operator L x on 'Hs by 

L x u :— u(x) — u(o). (2.2) 
Lemma 2.5. For any x 6 G°, one has \L x u\ < kS(u) 1 ^ 2 , where k depends only on x. 
Proof. Since G is connected, choose a path {*i}" =0 with xq = o, x n = x and c x . tX . , > for 
; = l,...,n. For k = {^ n i=l c~j x , J 1 , the Schwarz inequality yields 



\L x u\ = \u(x) - u(o)\ 



y / x " x - 1 ( u (xi) - u(Xi-t» 
tl V c Xi , Xi ., 



< k 2 &(u). 



Definition 2.6. Let v x be defined to be the unique element of 7Yg for which 

(v x , u)s = u(x) - u(o), for every u e < H? J . (2.3) 
This is justified by Lemma 2.5 and Riesz's lemma. The family of functions {v x } x ec° i s 
called the energy kernel because of Corollary 2.7. Note that v corresponds to a constant 
function, since (v„, u)g = for every u e "He. Therefore, this term may be ignored or 
omitted. 

Corollary 2.7. {vJ^go w a reproducing kernel for < Hg,- Thus, span{v x ) is dense in TYe. 

Proof. Choosing representatives with v x (o) = 0, it is trivial to check that (v x , v y )s = v x (v) = 
v v (x) and then apply Aronszajn's Theorem [Aro50]. The energy kernel is clearly total, so 
density is immediate. □ 

There is a rich modern literature dealing with reproducing kernels and their manifold 
application to both continuous analysis problems (see e.g., [ AD06, AL08, AAL08, B V03, 
Zha09]), and infinite discrete stochastic models. One novel aspect of the present work is 
the use of "relative" (to the reference vertex o) reproducing kernels. 

Remark 2.8. There is an alternative construction of 'Hg, via techniques of von Neumann 
and Schoenberg [vN32, Sch38b, Sch38a]. The natural notion of distance on G is the effec- 
tive resistance metric R, which is defined in terms of &. Because this metric is negative 
semidefinite, von Neumann's method gives an embedding cD : G — > Iri of the metric space 
(G,R) into a Hilbert space H in such a way that R(x,y) = ||d>(x) - C>(y)||^. In [JP09c, §6], 
it is shown that cD(x) i-» v x is a unitary isomorphism of Iri onto "Hs- 

Remark 2.9 (Probabilistic interpretation of v A ). The energy kernel {v v } is intimately related 
to effective resistance distance R(x, y). In fact, R(x, o) = v x (x) - v x (o) = S(v x ) and similarly, 
R(x,y) = S(v x - Vy). This is discussed in detail in [JP09b], but we give a brief summary 
here, to help the reader get a feeling for v A . For a random walk (RW) on a finite network 
started at the vertex y, let t x be the hitting time of x (i.e., the time at which the random 
walk first reaches x) and define the function 



u x (y) = ¥[t x < r 1 RW starts aty]. 

Here, the RW is governed by transition probabilities p(x,y) = c xy /c(x); cf. Remark 3.2. 
One can show that v x = R(x, d)u x is the representative of v x with v x (o) = 0. Since the range 
of u x is [0, 1], one has < v x (y) - v x (o) < v x (x) - v x (o) = R(x,o). Many other properties 
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of v x are similarly clear from this interpretation. For example, it is easy to compute v x 
completely on any tree. 

2.2. The finitely supported functions and the harmonic functions. 

Definition 2.10. For v e "7fg, one says that v has finite support iff there is a finite set 
F Q G° for which v(x) = k e C for all x t F. Equivalently, the set of functions of finite 
support in He, is 

span{<5. v ) = {u e dom£ '. u(x) = k for some k, for all but finitely many x e G°), (2.4) 

where 6 X is the Dirac mass at x, i.e., the element of He containing the characteristic func- 
tion of the singleton {x}. It is immediate from (1.14) that <5 A € Hs- Define Tin to be the 
closure of span{(5 A } with respect to £. 

Definition 2.11. The set of harmonic functions of finite energy is denoted 

Harm := {v € He, • Av(x) = 0, for all x e G°). (2.5) 
Note that this is independent of choice of representative for v in virtue of (1.10). 

Lemma 2.12. The Dirac masses {5 X } veC o form a reproducing kernel for A. That is, for any 
x e G , one has (6 X , w)g = Au(x). 

Proof. Compute (6 X , m)s = &{5 X , u) directly from formula (1.12). □ 

Remark 2.13. Note that one can take the definition of the Laplacian to be the operator A 
satisfying (S x , w)g = Au{x). This point of view is helpful, especially when distinguishing 
between identities in Hilbert space and pointwise equations. For example, if h e Harm, 
then Ah and the constant function 1 are identified in He because (u, A/t)g = (u, l)g = 0, 
for any u e *Hg. However, one should not consider a (pointwise) solution of Au(x) = 1 to 
be a harmonic function. 

Lemma 2.14. For any x e G°, Av x = 6 X - 5 . 

Proof. Using Lemma 2.12, Av x (y) = (d y , v v )g = 6 y (x) - 6 y (o) = (S x - 6 )(y). a 

Definition 2.15. A dipole is any v € satisfying the pointwise identity Av = S x — S y 
for some vertices x,y 6 G°. It is clear from Lemma 2.14 that the energy kernel consists of 
dipoles, and that v x - v y is always a finite-energy solution to Av = 5 X - 6 Y . 

The formula {5 x ,u)s = Au{x) of Lemma 2.12 is extremely important. Since Tin is 
the closure of span{(5 A ), it implies that the finitely supported functions and the harmonic 
functions are orthogonal. This result is sometimes called the "Royden Decomposition" in 
honour of the analogous theory established by Royden for Riemann surfaces, when it first 
appeared in [Yam79, Thm. 4.1]. However, the result is incorrect as stated there and the 
following corrected form may also be found in [Soa94, §VI] or [LP09, §9.3]. 

Theorem 2.16. 'Kg = TuKB'Harm. 

Proof. For all v 6 *Hg, Lemma 2.12 gives {6 X , v)g = Av(x). Since Tin = span{<5 v ), this 
equality shows v ± Tin whenever v is harmonic. Conversely, if (S x , v)g = for every x, 
then v must be harmonic. Recall that constants functions are in 'Hg. □ 

Corollary 2.17. span{(5 t ) is dense in iff 'Harm = 0. 
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Remark 2.18. Corollary 2.17 is immediate from Theorem 2.16, but we wish to empha- 
size the point, as it is not the usual case elsewhere in the literature and leads to unusual 
consequences, e.g., one may have 



More precisely, \\u - Yjxec t u(x)6 x \\b may not tend to as k — > oo, for any exhaustion {Gk\- 
Part of the importance of the energy kernel \v x ) arises from the fact that there is no other 
natural representing set for f Hg l , See also [Yam79, Lemma 3.1] for a similar interesting 
implication. 

Definition 2.19. Let f x = P^ m v x denote the image of v x under the projection to Tin. Simi- 
larly, let h x = P<HarmVx denote the image of v x under the projection to 'Harm. 

For future reference, we state the following immediate consequence of orthogonality. 

Lemma 2.20. With f x — P<f m v x , {f x } x ec° & a reproducing kernel for Tin, but f x J_ 'Harm. 
Similarly, with h x — P^ arm v x , {h x } xeC o is a reproducing kernel for T( arm, but h x ± Tin. 

While it may not be true that v y is in span{<5 A ) (or even in its closure), the following 
result shows that 5 y is always in span{v v ) when deg(>') < oo. 

Lemma 2.21. For any x e G°, 5 X — c(x)v x — Yiy~x c xy v y . 

Proof. Lemma 2. 12 implies (d x , u)g = (c(x)v x - Xy~x c xyV y , w)e for every u e 'Hs, so apply 
this to u — v z , z e G°. Since 8 X , v x e fig, it must also be that 2y~x c xy v y 6 "Hs- D 

Remark 2.22 (Real and complex-valued functions on G°). While we will need complex- 
valued functions for some results obtained via spectral theory, it will usually suffice to 
consider R-valued functions because of the following lemma. 

Lemma 2.23. The reproducing kernels v x , f x , h x are all H-valued functions. 

Proof. Computing directly, 



Then applying the reproducing kernel property, 

(v z , u)s = u(x) - u(o) = u(x) - u(o) = (v z , u)s- 
Thus (v^,u)s = (vz,u)e for every u 6 'Harm, and v z must be R-valued. The same compu- 



Definition 2.24. A sequence {m„} c 'He converges pointwise in 'Hs iff 3k € C such that 
u n (x) - u(x) — > k, for all x e G°. 

Lemma 2.25. If{u„} converges to u in S, then \u n ) converges to u pointwise in 'He- 
Proof Define w„ :- u„ - u so that ||w„||e — » 0. Then 





x,yeC° 



tation applies to f and h z . 



□ 



\W„(X) - W„(o)\ = \{V X ,W„) B \ < \\V X \\ S ■ ||w„||fi 

so that lim w„ exists pointwise and is a constant function. 



* 0, 



□ 
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3. The discrete Gauss-Green formula 



In Theorem 3.12, we establish a discrete version of the Gauss-Green formula which 
extends Proposition 1.9 to the case of infinite graphs; the scope of validity of this formula 
is given in terms of the space M of Definition 3. 1 . We are able to prove in Lemma 5.8 that 
the boundary term vanishes for elements of span{v A } and in Lemma 5.15 that it vanishes 
for finitely supported functions. Corollary 5.6 recovers the well-known fact that nontrivial 
harmonic functions cannot be in £ 2 (G°); see also the beginning of §4. 

A key difference between our development of the relationship between the Laplace 
operator A and the Dirichlet energy form S is that A is Hermitian but not necessarily self- 
adjoint in the present context. This is in sharp contrast to the literature on Dirichlet spaces 
in potential theory [Bre67,CC72] and the general theory of Dirichlet forms and probability 
[FOT94, BH91]. In fact, the "gap" between A and its self-adjoint extensions comprises an 
important part of the boundary theory for (G, c), and accounts for features of the boundary 
terms in the discrete Gauss-Green identity of Theorem 3.12. 

3.1. Monopoles and the domain of A. 

Definition 3.1. A monopole at x e G° is an element w x e Hs which satisfies Aw x (y) = 5 xy , 
where k e C and <5 lT is Kronecker's delta. When nonempty, the set of monopoles at the 
origin is closed and convex, so S attains a unique minimum here; let w always denote the 
unique energy-minimizing monopole at the origin. 

When He contains monopoles, let M x denote the vector space spanned by the monopoles 
at x. This implies that M x may contain harmonic functions; see Lemma 4.1. We indicate 
the distinguished monopoles 



where f x = Pj- m v x . (Corollary 4.4 below confirms that w v x = w f x for all x iff if Harm = 0.) 

Remark 3.2 (Monopoles and transience). The presence of monopoles in He is equivalent 
to the transience of the underlying network, that is, the transience of the simple random 
walk on the network with transition probabilities p(x,y) = c xy /c(x). To see this, note that 
if w x is a monopole, then the current induced by w x is a unit current flow to infinity with 
finite energy. It was proved in [ Yam79] that the network supports monopoles if and only if 
the Green kernel G exists (i.e., if and only if the random walk on the network is transient). 
This also appears in [Lyo83] in the form that the network is transient if and only if there 
exists a unit current flow to infinity; see [LP09, Thm. 2.10]. 

Let p„{x,y) denote the probability that the random walk started at x is aty after n steps. 
If y = (x = xq, x\, . . . , x„ = y) is a path from x to y, then define 



Let F(x, y) be the set of all finite paths from x to y, and let r„(x, y) be the subset of paths of 
length n. Then p\{x,y) — p(x,y) as defined just above, and in general, 



• x := v x + w„ and w f x :=f x + w, 



(3.1) 



n 




Pn(x,y)= J] P(r)- 



yeT„(x,y) 

Suppose the Green kernel G is defined in the usual manner: 
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0(x,y) = J^ Pn (x,y). (3.2) 

n=0 

Now let us denote the symmetrized Green kernel by g»(x,y) = @(x,y)/c(y), as in [Kig09a, 
§2]. Then g*(x,y) = g*(y, x) and there is a simple formula for g„ in terms of the (wired) 
resistance metric on G; see [Kig09a, (2.2)] and [JP09b]. For the moment, let us abuse 
notation and use w f x to denote the representative of w f x which vanishes at oo. Then the 
relationship between the monopoles and the Green kernel is given by 

w J x(y) = g*(x,y) = — — — . (3.3) 
c(y) 

Roughly speaking, the Green kernel is a reproducing kernel for the subspace of dom£ 
which consists of fi-limits of finitely supported functions. The monopole domain Ai con- 
tains an extension of the Green kernel which is a reproducing kernel for all of dom£. 

It is well-known that the transience/recurrence of the random walk is independent of the 
vertex at which the walk is started; equivalently, Q(x, y) is finite for some x and y if and 
only if it is finite for all x, y. The corresponding statement for monopoles is also clear from 
(3.1): there exists a monopole at some vertex o if and only if there exists a monopole at 
any other vertex x. 

Remark 3.3. Note that w a e Tin, whenever it is present in "Kg, and similarly that w{ x is the 
energy-minimizing element of Ai x . To see this, suppose w x is any monopole at x. Since 
w x 6 'Hs, write w x = f + h by Theorem 2.16, and get &{w x ) = £(/) + S(lt). Projecting 
away the harmonic component will not affect the monopole property, so w x = P<r m w x is 
the unique monopole of minimal energy. Also, w corresponds to the projection of 1 to Do; 
see §4.1. 

Definition 3.4. The dense subspace of 94$ spanned by monopoles (or dipoles) is 

Ai := span) \> x } xeC ° + spanjw^, w^} ieG o . (3.4) 
Let Am be the closure of the Laplacian when taken to have the dense domain At. 

Note that Ai = span{v c ) when there are no monopoles (i.e., when all solutions of of 
Aw = 6 X have infinite energy), and that Ai = span{w^., w f x ) when there are monopoles; see 
Lemma 4. 1 . 

The space Ai is introduced as a dense domain for A (see Remark 3.8) and for its use 
as a hypothesis in our main result, that is, as the largest domain of validity for the discrete 
Gauss-Green identity of Theorem 3.12. Note that while a general monopole need not be in 
dom Am (see [JP09c, Ex. 13.8 or Ex. 14.39]), we show in Lemma 3.7 that it is always the 
case that it lies in dom A^. 

Definition 3.5. A Hermitian operator S on a Hilbert space H is called semibounded iff 

(v, S v) > 0, for every v € D, (3.5) 
so that its spectrum lies in some halfline [k, oo) and its defect indices agree. 

Lemma 3.6. Am is Hermitian; a fortiori, Am is semibounded. 
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Proof. Suppose we have two finite sums u — £ a x w x and v = £ b y w y , writing w x for w v x or 
w f x . We may assume that o appears neither in the sum u nor for v; see Definition 2.6. Then 
Lemma 2.12 gives 

<m,Av) £ = ^ ~a x b y (w x , Aw y )e = ^a7fo v (w v ,£ v > £ = ^a^Aw^Cy) = ^ajfo,.^,,. 

Of course, (Am, v)g = 2 ~a~xby5 xx exactly the same way. The argument for linear combina- 
tions from \v x ] is similar, so Am is Hermitian. Then 

(u, Au) s = ^ a~ K a y 5 xy = ^ |a t | 2 > 

x,y x 

shows Am is semibounded. The argument for {v x } is similar. □ 

Since A agrees with Am pointwise, we may suppress reference to the domain for ease 
of notation. When given a pointwise identity Au = v, there is an associated identity in 
but the Lemma 3.7 shows that one must use the adjoint. 

Lemma 3.7. For u, v e Tig, Au — v pointwise if and only if v — A M u in TYg. 

Proof. We show that u e domAj^ for simplicity, so let tp e span{v A ) be given by tp = 
YIi=\ fliVv,; the proof for tp e span{w^, w f x ] is similar. Then Lemma 2.12 and Lemma 2.14 
give 

n n 

{Aip, u)e = ^ a.i(6 Xi - 6 , u) E = ^ a^AuiXi) - Au(o)). 

!=1 1=1 

Since Am(x) = v(x) by hypothesis, this may be continued as 

n n 

(A(p, u) s = 2] «f(v(JCj) - v(o)) = ^ a,(vj,, v) 6 = {ip, v> £ . 

!=1 !=1 

Then the Schwarz inequality gives the estimate \{A<p,u)&\ = \{(p,v)s\ < \\(p\\s\M\s, which 
means u € domA^. The converse is trivial. □ 

Remark 3.8 (Monopoles give a reproducing kernel for ran Am). Lemma 3.7 means that 

{w x , Au)s = (S x , u)e, for all u e dom A^- (3.6) 
for every w x e M x . Combined with Lemma 2.12, this immediately gives 

(w x ,Au) e = Au(x). (3.7) 
If {wxLgg i s a collection of monopoles which includes one element from each At*, then 
this collection is a reproducing kernel for ran A;m. Note that the expression Au(x) is defined 
in terms of differences, so the right-hand side is well-defined even without reference to 
another vertex, i.e., independent of any choice of representative. 

In Definition 3.4, we give a domain M for A which ensures that ranAyn contains all 
finitely supported functions and is thus dense in Tin. However, even when A is defined 
so as to be a closed operator, one may not have Tin C ran A; in general, the containment 
ran(5 cl °) C (ran S ) cl ° may be strict. The operator closure S cl ° is done with respect to the 
graph norm, and the closure of the range is done with respect to S. We note that [MYY94, 
(G. 1)] claims that the Green function is a reproducing kernel for all of Tin. In our context, 
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at least, the Green function is a reproducing kernel only for ran A, where A has been chosen 
with a suitable dense domain. In general, the containment ran A c Tin may be strict. In 
fact, it is true that ran A^ c Tin, and even this containment may be strict. Note that w{ x is 
the only element of M x which lies in (ran Am) c1 °, and it may not lie in ran Am- 

A different choice of domain for A can exacerbate the discrepancy between ran A and 
Tin: if one were to define Ay to be the closure of A when taken to have dense domain 
"V := span{v A ), then ran Ay is dense in Tinn, the £-closure of spanj^ - 6,,}. However, it 
can happen that Titij is a proper orthogonal subspace of Tin (the fi-closure of span{5 A }). 
This is discussed further in [JP09c, §8]; an example of / € Tin\ := Tin Q Tin2 is com- 
puted in [JP09c, Ex. 14.35]. The domain of A can thus induce a refinement of the Royden 
decomposition: 

"Hs = Tin \ © Tiriz © "Harm. 

See Theorem 2.16 and the comment preceding it. This highlights the importance of the 
choice of M as the domain for A. 

3.2. Relating A to £. 

Definition 3.9. If H is a subgraph of G, then the boundary of H is 

bdH :={xeH 3yeH c ,y ~ x). (3.8) 
The interior of a subgraph H consists of the vertices in H whose neighbours also lie in H: 

intff := {x € H \y ~ x => y e H] = H \ bdH. (3.9) 
For vertices in the boundary of a subgraph, the normal derivative of v is 

2 c ^ v{x) ~ V W)> for xebdH. (3.10) 

yefl 

Thus, the normal derivative of v is computed like Av(jt), except that the sum extends only 
over the neighbours of x which lie in H. 

Definition 3.9 will be used primarily for subgraphs that form an exhaustion of G, in the 
sense of Definition 1.5: an increasing sequence of finite and connected subgraphs [Gk], 
so that Gk c Gk+i and G = {JGu- Also, recall that 2bdc : - limfc-xw ZbdGi from Defini- 
tion 3.10. 

Definition 3.10. A boundary sum is computed in terms of an exhaustion [Gk] by 

bdC bdG t 

whenever the limit is independent of the choice of exhaustion, as in Definition 1.6. 

Remark 3.11 (bdG as a measure space). The boundary bdG is formalized as a measure 
space in [JP09a, JP09d]; see also [JP09c, §7]. Just as the Martin boundary is a support 
set for a representing measure for the nonnegative harmonic functions, and the Poisson 
boundary is a support set for a representing measure for the bounded harmonic functions, 
it is shown in [JP09a] that bd G is a support set for a representing measure for the harmonic 
functions of finite energy. For more about Martin and Poisson boundaries, see [Saw97, 
WoeOO] and the references therein. 
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1 Theorem 3.12 (Discrete Gauss-Green Formula). Ifu e and v e M, then 

< M ,v) £ = 2" Av + Z"^- (3 ' 12) 

G° bdG 

2 Proof. It suffices to work with R-valued functions and then complexify afterwards. By the 

3 same computation as in Proposition 1.9, we have 

i 2 c xy (u(x)-u(y))(v(x)-v(y))= J] u(x)Av(x) + J] ( 3 - 13 ) 

x, veCi teint (.' .tebd (.'. 

4 Taking limits of both sides as k — > oo gives (3.12). It remains to see that one of the sums 

5 on the right-hand side is finite (and hence that both are). Note that if v = w z is a monopole, 

6 then 

^ u(x)Av(x) = ^ u(x)6 z (x) = u(z). 

xeC° .veC° 

7 This is obviously independent of exhaustion, and immediately extends to v e M. □ 

8 The following boundary representation of harmonic functions is the focus of [JP09a]. 

9 Corollary 3.13 (Boundary representation of harmonic functions). For all u e "Harm, 

u(x) = ^wf^ + u(o). (3.14) 

bdC 

10 Proof. Lemma 2.23 and (2.3) imply u(x) - u(o) = (v x , m)s = (u, v x )s = ZbdG M ffr ■ D 

11 Lemma 3.14. For all u e domAyn, 2c° Am = — 2bdc f^- Thus, the discrete Gauss-Green 

12 formula (3.12) is independent of choice of representatives. 

13 Proof. On each Gk, each edge appears twice in the sum (once with each sign/orientation) 

14 and so 

J] A«(jc) + J] J£(x) = ^ c A ,,( M (x) - «(y)) = 0. 

.reint G s .tebd C ( x,y£G t 

15 To check a different representative, use the first part to compute 




□ 



16 Remark 3.15. It is clear that (3.12) remains true much more generally than under the 

17 specified conditions; certainly the formula holds whenever 2.vec° l«(x)Av(x)| < oo. Un- 

18 fortunately, given any hypotheses more specific than this, the limitless variety of infinite 

19 networks almost always allow one to construct a counterexample; i.e. one cannot give a 

20 condition for which the formula is true for all u e < H&, for all networks. To see this, suppose 

21 that v = Y^\L\ a (w x , with each w Xl a monopole at the vertex x,. Then 

oo 

^ u(x)Av(x) = ^ a.iu(xi), 

A-eG° i=l 
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and one would need to provide a condition on sequences {a,} that would ensure a,M(jc,) 
is absolutely convergent for all u e 'Hg. Such a hypothesis is not likely to be useful (if it 
is even possible to construct) and would depend heavily on the network under investiga- 
tion. Nonetheless, the formula remains true and even useful in many specific and general 
contexts. For example, it is clearly valid whenever v is a dipole, including all those in 
the energy kernel. We will also see that it holds for the projections of v x to Tin and to 
Harm. Consequently, for v which are limits of elements in M, we can use this result in 
combination with ad hoc arguments. 

Remark 3.16. We refer to 2bdc as me "boundary term" by analogy with classical PDE 
theory. This terminology should not be confused with the notion of boundary that arises in 
the discussion of the discrete Dirichlet problem, where the boundary is a prescribed set of 
G°. As the boundary term may be difficult to contend with, it is extremely useful to know 
when it vanishes, for example: 

(i) when the network is recurrent (Lemma 4.6), 

(ii) when v is an element of the energy kernel (Lemma 5.8), 

(iii) when u, v, Am, Av lie in ( 2 (Lemma 5.12), and 

(iv) when either u or v has finite support (Lemma 5.15). 

In fact, Lemma 4.6 expresses the fact that it is precisely the presence of monopoles that 
prevents the boundary term from vanishing. 

4. More about monopoles and the space M 

This section studies the role of the monopoles with regard to the boundary term of 
Theorem 3.12, and provides several characterizations of transience of the network, in terms 
the operator- theoretic properties of Am- 

Note that if h e "Harm satisfies the hypotheses of Theorem 3.12, then£(/i) = Zbdc^f^- 
In Theorem 4.6 we show that &(u) = £ G o uAu for all u e "Hg iff the network is recurrent. 
With respect to "Kg = Tin © fiarm, this shows that the energy of finitely supported func- 
tions comes from the sum over G°, and the energy of harmonic functions comes from the 
boundary sum. However, for a monopole w x , the representative specified by w x {x) = 
satisfies &(w) = ZbdG w |^ but the representative specified by w x (x) = S(w x ) satisfies 
fi(vv) = 2g° wAw. Roughly speaking, a monopole is therefore "half of a harmonic func- 
tion" or halfway to being a harmonic function. A further justification for this comment is 
given by Corollary 4.4: the proof shows that a harmonic function can be constructed from 
two monopoles at the same vertex. A different perspective one the same theme is given in 
Remark 4. 14. The general theme of this section is the ability of monopoles to "bridge" the 
finite and the harmonic. 

Lemma 4.1. When the network is transient, M contains the spaces span{v A }, span{/ x }, and 
span{/i v }, where f x = P Ti „v x and h x = P, Han „v x . 

Proof. The first two are obvious, since v x = w\ - w and f x = w f x — w by Definition 3.1. 
For the harmonics, note that these same identities give 

w v x - w — v x = f x + h x = w{ x - w + h x , 
which implies that h x = w\ — w f x . (Of course, w x = w f x when fiarm = 0.) □ 

Theorem 4.2 ( [Soa94, Thm. 1.33]). Let u be a nonnegative function on a recurrent net- 
work. Then u is superharmonic if and only if u is constant. 
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Corollary 4.3. If Harm + 0, then there is a monopole in "He,. 

Proof. If h 6 'Harm and h + 0, then h = hi - I12 with /?, e Harm and h, > 0, by [Soa94, 
Thm. 3.72]. (Here, hi > means that hi is bounded below, and so we can choose a 
representative which is nonnegative.) Since the /?, cannot both be 0, Theorem 4.2 implies 
the network is transient. Then by [Lyo83, Thm. 1], the network supports a monopole. □ 

Corollary 4.4. Harm + iff there is more than one monopole at one (i.e. every) vertex x. 

Proof. If He contains a monopole w x + w v x , then h := w v x - w x is a nonzero element of 
Harm. Conversely, if Harm + 0, then there is a monopole w e He by Corollary 4.3. For 
a nonzero element h e Harm, w + h is also a monopole. □ 

Definition 4.5. The phrase "the boundary term is nonvanishing" indicates that (3.12) holds 
with nonzero boundary sum when applied to {u, v)e, for every representative of u except 
one; namely, the one specified by u(x) = (u, w v x )e- 

Recall from Remark 3.2 that the network is transient iff there are monopoles in He- 

Theorem 4.6. The network is transient if and only if the boundary term is nonvanishing. 

Proof. (=>) If the network is transient, then as explained in Remark 3.2, there isawe He 
with Aw = S z . Now let w z := iV,„w so that for any u e dom Am, (3.12) 



{u,w z )e = u(z) + 2^ u Tfa- 



bdG 

It is immediate that 2bdG = if and only if the computation is done with the repre- 
sentative of u specified by u{z) = (u, w z )e- 

(<=) Suppose that there does not exist w e He with Aw = S z , for any z e G°. Then 
Ai = span{v x } as discussed in Definition 3.1. Therefore, it suffices to show that 

(u, v x )e = 2] uAv x , 

but this is clear because both sides are equal to u(x) - u(o) by (2.3) and Lemma 2.14. □ 

The next theorem shows that monopoles may be constructed explicitly as weak-* limits 
in the Hilbert space He- Nonuniqueness of weak-* limits may lead to nonuniqueness of 
monopoles at a given point x (provided Harm + 0; see Corollary 4.4). One would like 
to argue as follows (as suggested by a referee): for s > 0, let D s be the multiplication 
operator corresponding to diagonal matrix with entries < 1 . Then with p(x, y) = 

Cxy/cix) giving transition probabilities as before, and defining the heat operator pointwise 
by Pu(x) = ^yec° P( x > y) u (y)> one can see tnat tne row sums of the corresponding matrix 
satisfy 



{D s Pl)(x) = V C(X) = C(X) < 1 , for each fixed x e G 

*— l n e + c(x) c(x) e + c(x) 

veG° 

Consequently, 







- . CO , DO 

(e + A)- 1 = (I - D £ P) 1 = V(D B F) n -^U - V p(x,y) n . (4.1) 

S + C E + C *—i C 
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It is well known that the network is transient if and only if the entries of the matrix 
2j«^o P( x <y)" corresponding to the Green operator are finite. However, this runs into diffi- 
culties in the case when c is unbounded, and so we make use of the spectral theorem. 

Theorem 4.7. The network is transient if and only if f k := (e k + Am) \ 6 x is weak-* con- 
vergent for some sequence s k — > 0. 

Proof. We show that there is a monopole if and only if there is sequence \s k } with s k — > 
and sup^ \\(s k + AvO^Hs < B < oo. 

(=>) Let Am be any self-adjoint extension of Am whose spectrum lies in [0, oo), and let 
E(dA) be the corresponding projection-valued measure. For concreteness, one may take 
the Friedrichs extension, but this is not necessary; A commutes with conjugation, and so a 
theorem of von Neumann implies that such an extension exists. Then 



E(dA)u, (4.2) 

o e + A 

where we use the notation R e := (s + Am)" 1 for the resolvent. Note that AmR e £ (Am/? s )* = 
R* S A* M = R s Am- On the other hand, Am £ A^ and therefore R b Am £ B e A^. Combining 
these gives AmR e £ B. e A' M . Now we apply this and (4.2) to u — A*w to get 

r°° a 

f k = (s k + Aa,) -1 ^- = (£k + Am)" 1 Av,w = A M (s k + A M T l w = -E(dA)w. 

Jo e k + A 

Note that R E is bounded, and so w e domR e automatically. This integral implies 



II/* 



\l< j ) \\E(dA)w\\l< \ \\E(dA)w\\ 2 s = |M||. 

Jo \Ek + AJ Jo 



Thus we have supj. ||(s* + A) l 8 x \\ e = sup ||/*|| e < B - \\w\\ E < oo. 

(<=) We show the existence of a monopole at x. Since e k f k + Af k = 6 X , the bound 
SU P WfkWs ^ B implies that 

IIAf* - S x \\ s = \\s k f k \\ < s k B 0. 
Let w be a weak-* limit of \f k }. Then for if e dom Am, 

(A<p,w)s = ]im(Aip,f k ) s = lim(ip, Af k ) s = lim<^,<5 r - s k f k ) £ = {ip,8 x ) E , 

k — »0O k — JOO k — >00 

so that w is a monopole at x. □ 

Lemma 4.8. On any network, (ran Am) c '° c Tin and hence 'Harm c ker A^ ( . 

Proof. If v e M, then clearly Amv e Tin. To close the operator, we consider sequences 
{u n } c M. which are Cauchy in £, and for which {Am,,} is also Cauchy in £, and then include 
u := limM,, in dom Am by defining Amm := lim Am«„. Since /„ := Am m « has finite support 
for each «, the £-limit of {/„} must lie in Tin. Since Tin is closed, the first claim follows. 
The second claim follows upon taking orthogonal complements. □ 

Theorem 4.9. The network is transient if and only ;/(ran A^) rf = Tin. 
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1 Proof. (=>) If the network is transient, we have a monopole at every vertex; see Re- 

2 mark 3.2. Then any u e span{<5 x } is in ranAj^ because the monopole w x is in domAvi, 

3 and so Tin c ran A^. The other inclusion is Lemma 4.8. 

4 (<=) If S x € ran Aw for some x e G°, then A^vv = S x for w e domAvf c domfi and so 

5 w is a monopole. Then the induced current dw is a unit flow to infinity, and the network is 

6 transient, again by [Lyo83]. □ 

7 4. 1 . Comparison with the grounded energy space. There are some subtleties in the re- 
a lationship between «H fi and D as discussed in [LP09] and [KY89, KY84, MYY94, Yam79], 

9 so we take a moment to give details. We have attempted to match the notation of these 

10 sources. 

11 Definition 4.10. The inner product 

(u, v) := u{o)v(o) + (u, v) s . 

12 makes domfi into a Hilbert space D which we call the grounded energy space. Let Do be 

13 the closure of span{£ A } in D and let HD be the space of harmonic functions in D. 

u Throughout this section (only), we use the notation u :- u(o), for u e D. 

15 Definition 4.11. With regard to D, we define the vector subspace 

M~ := [u € D • Am = -u 5 ). (4.3) 

16 Note that M~ contains the harmonic subspace 

HD C := {u € D • Am = and u B = 0). (4.4) 

17 The previous definition is motivated by the following lemma, 
is Lemma 4.12. D^ = M~ and hence D = Do © M~. 

19 Proof. With m := m(o), we have u 6 Dq iff u ± span{5 v ), which means that 

= <m, 8 X ) = u 5 x {o) + <m, 8 x )s = u 5 xo + Au(x), Vx € G°, (4.5) 

20 which means Am = —u 8 . □ 

21 Let us denote the projection of D to Do by Pd and the projection to Dq by P^. 

22 Remark 4. 1 3. The constant function 1 decomposes into a linear combination of two monopoles: 

23 let v = Pd 1 and u = ^ 1 = 1 - v, and observe that Am = -u 5 by Lemma 4.12 and that 

24 Av = A(l - m) = -Am = u 5 , so u„ — 1 - v a gives Av = (1 - v )5 . In general, the constant 

25 function kl decomposes into v = Pi) kl and u = kl, where 

Av = (k — u )5 and Am = —u 8 . 

26 With respect to the decomposition D = Do © M~, given by Lemma 4.12, there are two 

27 monopoles w„ e Do and 6 M~ (which may be equal) such that 1 = u w ( ^ - u w . 

28 When one passes from D to 'Kg by modding out constants, these components of 1 add 

29 together to form (possibly constant) harmonic functions. An example of this is given in 

30 Example 6.8. 
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Consequently, Lemma 4.12 yields a short proof of [LP09, Exc. 9.6c]: Prove that the 
network is recurrent iff 1 e Do. To see this, observe that if u is the projection of 1 to 
Dq , then u + iff there is a monopole. This result first appeared (in more general form) 
in [Yam77, Thm. 3.2]. 

Remark 4.14. Despite the fact that Theorem 2.16 gives "Hg = Tin © Harm, note that 
D + Do © HD. This is a bit surprising, since "Hg = D/Cl, etc., and this mistake has been 
made in the literature, e.g. [Yam79, Thm. 4.1]. The discrepancy results from the way that 
1 behaves with respect to Pd ; this is easiest to see by considering 

D +*:={/ + *l-/e D , ifceC}, k + 0. 
If the network is transient and / e Do + k, k + 0, then / = g + kl for some g e Do, and 

f = (g + kP Do l) + kP^l 

shows / i Do. Nonetheless, it is easy to check that Do + k is equal to the o-closure of 
span 6 X + k, and hence that (D + C1)/C1 = Tin. This appears in [LP09, Exc. 9.6b]. 
Similarly, note that for a general h e HD, one has h = P^Ji + kl, so that h i Dq. 

We conclude with a curious lemma that can greatly simplify the computation of monopoles 
of the form Pd 1; it is used in Example 6.8. In the next lemma, u = u(o), as above. 

Lemma 4.15. Let u e Dq. Then u = P^l if and only if u a = S(u) + u 2 D € [0, 1). 

Proof. From ||m||^ + ||1 - uW 2 , — \\\\\^ - 1, one obtains S(u) - u a + \u \ 2 = 0. From (u, 1 - 
u)„ = 0, one obtains &(u) -u^+ \u \ 2 = 0. Combining the equations gives u a = u„ = 
i(l ± yi - 4£(u)), so that u e [0,1]. However, u + 1 or else &(u) = would imply 
1 € Dy in contradiction to (4.5). The converse is clear. □ 

Remark 4.16. The significance of the parameter u„ is not clear. However, it appears to 
be related to the overall "strength" of the conductance of the network; we will see in 
Example 6.8 that u„ « 1 corresponds to rapid growth of c near oo. Also, it follows from the 
Remark 4. 1 3 and Lemma 4.15 that u = corresponds to the recurrence. There is probably 
a good interpretation of u in terms of probability and/or the speed of the random walk, but 
we have not yet determined it. The existence of conductances attaining maximal energy 
&{P^\) = | is similarly intriguing, and even more mysterious. Example 6.8 shows that 
the maximum is attained on (Z, c") for c — 2. 

5. Applications and extensions 

In §5.1, we use the techniques developed above to obtain new and succinct proofs of 
four known results, and in §5.2 we give some useful special cases of our main result, 
Theorem 3.12. 

Definition 5.1. For an infinite graph G, we say u(x) vanishes at oo iff for any exhaustion 
{Gk}, one can always find k and a constant C such that \\u{x) - C\\oo < £ for all x $ Gk- 
One can always choose the representative of u 6 Hs so that C — 0, but this may not be 
compatible with the choice u(o) = 0. 

Definition 5.2. Say y = (xq, x\,xz, . . .) is a path to oo iff % i ~ for each i, and for any 
exhaustion {Ga) of G, 
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Vk, 3N such that n>N=* 



x„ ( G k . 



(5.1) 



5.1. More about Tin and 'Harm. The next two results are almost converse to each other, 
although the exact converse of Lemma 5.3 is false; see [JP09c, Fig. lOorEx. 14.16]. How- 
ever, the converse does hold almost everywhere with respect to the usual (Kolmogorov) 
measure on the space of paths starting at o, by [ ALP99, Cor. 1 .2] (see [Car72] for con- 
struction of the measure). Lemma 5.3 is also related to [Soa94, Thm. 3.86], in which the 
result is stated as holding almost everywhere with respect to the notion of extremal length. 

Lemma 5.3. If u e He and u vanishes at oo, then u e Tin. 

Proof. Let u — f + h e Us vanish at oo. This implies that for any exhaustion {Ga ) and any 
s > 0, there is a A: and C for which \\h{x) - CIU < s outside Gk. A harmonic function can 
only obtain its maximum on the boundary, unless it is constant, so in particular, s bounds 
\\h(x) - C\\oo on all of Gk- Letting s — > 0, h tends to a constant function and u — f. □ 

Lemma 5.4. If h e Harm is nonconstant, then from any xq e G°, there is a path to infinity 
y — (xq, x\ , . . . ), with h(xj) < h(xj+i)for all j = 0, 1, 2, . . . . 

Proof. Abusing notation, let h be any representative of h. Since h(x) = 2y~.v iTfj^OO - 
sup y ^ v h(y) and h is nonconstant, we can always find y ~ x for which h(y{) > /i(xo). This 
follows from the maximal principle for harmonic functions; cf. [LP09, §2.1], [LPW08, 
Ex. 1.12], or [Soa94, Thm. 1.35]. Thus, one can inductively construct a sequence which 
defines the desired path y. Note that y is infinite, so the condition h(xj) < h(xj+i) eventually 
forces it to leave any finite subset of G°, so Definition 5.2 is satisfied. □ 

It is instructive to prove the contrapositive of Lemma 5.3 directly: 
Lemma 5.5. If h e Ttarm \ {0}, then h has at least two different limiting values at oo. 

Proof. Choose x € G° for which h x = P<HarmV x £ 'Kb is nonconstant. Then Lemma 5.4 
gives a path to infinity y\ along which h x is strictly increasing. Since the reasoning of 
Lemma 5.4 works just as well with the inequalities reversed, we also get 72 to 00 along 
which h x is strictly decreasing. This gives two different limiting values of h x , and hence h x 
cannot vanish at 00. □ 

Corollary 5.6. Ifh e Warm is nonconstant, then h t £ p (G°)for any 1 < p < 00, 

Proof. Lemma 5.5 shows that no matter what representative is chosen for h, the sum \\h\\ p = 
2 veg° \K x )\ p nas the lower bound 2 x£ f e '' - S ''\F\, for some infinite set F c G°. □ 

5.2. Special cases of the Discrete Gauss-Green formula. In this subsection, we establish 
that the boundary term of Theorem 3. 12 vanishes for vectors in span{v A ] in Lemma 5.8, and 
that A is Hermitian when its domain is correctly chosen, in Corollary 5.10. These results 
continue the theme of Theorem 4.6. 

Lemma 5.7. For u € span{v x ], 2.vec° & u ( x ) - 0. 

Proof. For a finite sum u = 2 a y v y , the result follows by interchanging finite sums: 



2^ Av y (x) = 2^ a y (5 y - 6 )(x) = 2^ a.v - 2^ 



a x = 0. 



□ 



xeC° x 
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Physically, £ Am is the net divergence of the current passing through the network. Thus, 
Lemma 5.7 can be rephrased as saying that elements of span{v x ) are "balanced"; compare 
to [Soa94, p. 45]. Lemma 5.7 is false for u = w x e M x and may also fail for u in the 
closure of span{v v ) (with respect to S or the graph norm of A). 

Lemma 5.8. If v e span{v v ), then (u, v)g = £ reG o u(x)Av(x). 

Proof. It suffices to consider v = v x , whence 

^ u(y)Av x (y) = ^ u(y)(S x - S )(y) = u(x) - u(o) = {u, v x ) E , 
c° G° 

by Lemma 2.14 and the reproducing property of Corollary 2.7. □ 

Note that the formula in Lemma 5.8 may look odd because the right-hand side appears 
to depend on a choice of representatives, but this not the case by Lemma 5.7. However, 
Lemma 5.8 is false for v e M \ span{v t ). 

When deg(.Y) < oo for all x e G°, the next result follows immediately from Lemma 5.8 
and Lemma 2.21. However, it is true even when there are vertices of infinite degree. 

Theorem 5.9. For u, v e span{v A ), 

<m,A\'} £ = V Au(x)Av(x). 

jceG" 



(5.2) 



Proof. Let u € span{v A ) be given by the finite sum u = Yi X %xV x - We may assume the sum 
does not include o (see Definition 2.6). Then 



(5.3) 



Now we have 



(U, Au) E = Y £x£y(Vx, Av y > 6 = Y ZxtyiVx, 6 y - 6„) E . 
x,y x,y 

Since it is easy to compute (v x , 6 y - 5 )s = 6 xy + 1 (Kronecker's delta), we have 



(u, Au) s = J] te y {6 v + 1) = 2 \tA 2 + 2 & = J] |A«MI 2 + 2 Am(x) 

A',V X X X X 

by (5.3). The second sum vanishes by Lemma 5.7, and (5.2) follows by polarizing. 
Corollary 5.10. Am satisfies 



0<Y |A«(x)| 2 < (u, Au) s 



< oo. 



xeC° 



(5.4) 



(5.5) 



Proof. For u e At, we can repeat the proof of Theorem 5.9 up to (5.4) to obtain the central 
inequality of (5.5); the inequality appears because the second sum may not vanish in this 
case. Now let u € domAyn and choose {«„} c M. with Hindoo \\u n - - lim,,^^ \\Au„ - 
Au\\s = 0. Then the result follows by Fatou's lemma as in [Mal95, Thm. 1.7.7], □ 
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Remark 5.11. The notation u e t x means J^xec" \ u ( x )\ < 00 an d the notation u e 1 2 means 
Tjxec" \ u ( x )\ 2 < 00 • When discussing an element u of 7Ye, we say u lies in ( 2 if it has a 
representative which does, i.e., if u + k e { 2 for some k e C. This constant is clearly 
necessarily unique on an infinite network, if it exists. 

The next result extends Proposition 1.9 and is a partial converse to Theorem 3.12. All 
that is required for the computation in the proof of Lemma 5. 12 is that uAv e l x , which is 
certainly implied by u, Av e I 1 . However, this would not suffice to show u, v e dom£. 

Lemma 5.12. Ifu, v, Am, Av e I 1 , then (u, v)g = 2. ve g° u(x)Av(x), and u,v € domfi. 

Proof. If u, Av e I 1 , then mAv € t , and the proof of Proposition 1.9 is still valid: the 
absolute convergence of Yjxec u(x)Av(x) justifies the rearrangement. Substituting u for v 
in this formula gives u e dom£, and similarly for v. □ 

Lemma 5.13. Let w z e Ti^ be a monopole at z. Then w z e ran(I + A^) if and only if there 
is a function u e satisfying u + Am = on G° ' \ {z}. 

Proof. (=>) Since w z = v + Av for some v e Ws, set u :- v - w z . Then for x ± z, it is easy 
to check m(x) + Am(x) = 0. (<=) Set v :- w z + au for a := -1/(m(z) + Am(z)). Then it is easy 
to check that v(x) + Av(x) = w z (x). a 

If there are vertices of infinite degree in the network, then it does not necessary follow 
that span{£ r } c span{v A ), or that span{<5 A ) c M. However, we do have the following result. 

Definition 5.14. Let T := span{(5 A } ie G° denote the vector space of finite linear combina- 
tions of Dirac masses, and let Af be the closure of the Laplacian when taken to have the 
domain T ■ 

Note that T is a dense domain only when Harm = 0, by Corollary 2.17. Again, since A 
agrees with Af pointwise, we may suppress reference to the domain for ease of notation. 
The next result extends Proposition 1.9 to networks with vertices of infinite degree. 

Theorem 5.15. Ifuorv lies in dom Af, then (u, v)g = Yixec u(x)Av(x). 

Proof. First, suppose u e domA^ and choose a sequence {u„} C span{i5. v } with ||w n -w||s — > 
0. From Lemma 2.12, one has (6 x ,v)s = Av(x), and hence (u n ,v)e = 2.vec° u n (x)Av(x) 
holds for each n. Define M := sup{||«„||e}, and note that M < 00, since this sequence is 
convergent (to ||u||g). Moreover, \{u n , v)g| < M ■ ||v||g by the Schwarz inequality. Since u n 
converges pointwise to u in "Kg by Lemma 2.25, this bound will allow us to apply Fatou's 
Lemma (as stated in [Mal95, Lemma 7.7], for example), as follows: 



where we have used the hypothesis followed by u„ e span{i5 v ). Note that the sum over G 
is absolutely convergent, as required by Definition 1.5. 

Now suppose that v € dom Ayr and observe that this implies v e Tin also. By Theo- 
rem 2.16, one can decompose u = f + h where / = P Tm u and h = P. Hcmn u, and then 




v) s = (/, v>£ + (h, v) B = (f, v)e, 
since h is orthogonal to v. Now apply the previous argument to (/, v)g. □ 
For networks of finite degree, Theorem 5 . 1 5 follows from Theorem 3 . 1 2 by Lemma 2.21. 
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6. Examples 

In this section, we introduce the most basic family of examples that illustrate our tech- 
nical results and exhibit the properties (and support the types of functions) that we have 
discussed above. After presenting some basic examples, we prove some theorems regard- 
ing the properties of these examples. 

Networks similar to Example 6. 1 have been discussed elsewhere in the literature (for 
example, [KY89, Ex. 3.12, Ex. 4.9] and [KY84, Ex. 3.1, Ex. 3.2]), but the authors appear 
to assume that A is self-adjoint. This is not generally the case when c is unbounded; 
in fact, the Laplacian is not self-adjoint for Example 6.3 or Example 6.2. The proof is 
unfortunately beyond the scope of this paper; see [JP09c, Ex. 14.36 and Ex. 14.39] for 
further discussion and the explicit computation of defect vectors. 

Example 6.1 (Integer networks). Let (Z, c) denote the network with integers for vertices, 
and with conductances defined by c. We fix o — 0. 

These networks are more interesting when c grows fast enough to ensure, for example, 
that 2 c X y < 00 ■ I n this case, it is helpful to keep the following more concrete model in 
mind, especially if one hopes for tractable computations. 

Example 6.2 (Geometric integer model). For a fixed constant c > 1, let (Z, c") denote the 
network with integers for vertices, and with geometrically increasing conductances defined 
by c„_i,„ = c max lM>l«-i|) so m at the network under consideration is 



Again, we fix o — 0. Theorem 6.4 shows that Harm + for this network, and Lemma 6.5 
exhibits dipoles, monopoles, and a harmonic function on Example 6.2. 

Example 6.3 (Geometric half-integer model). It is also interesting to consider (Z + , c"), as 
this network supports a monopole, but has 'Harm = 0. 



Lemma 6.5 exhibits dipoles and a monopole for this example, but this network does not 
support harmonic functions. 

For k = 2,3,..., the network (Z+,k") can be thought of as the "projection" of the 
homogeneous tree of degree k (Tjt, il) under a map which sends x ton eZ iff there are n 
edges between x and o. 

Theorem 6.4. Harm + for (Z, c) iff 2 c ty < 00 • ^ n m ' s case > Harm is spanned by a 
single bounded function. 

Proof. (=s>) Fix m(0) = 0, define u{\) - ^- and let u(n) be such that 

1 

u(n) — u(n — 1) = , Vn. (6.1) 

Now u is harmonic: 



Au(n) = c„_i „(M(n) - u(n - 1)) - c n n+ \(u{n + 1) - u{n)) 
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and u is of finite energy 

S(u) = V c B _i fB (n(n) - u(n - l)) 2 = Y — !— < oo. 

^ c„_i,„ 

neZ «eZ 

Note that once the values of m(0) and m(1) are fixed, all the other values of u(n) are deter- 
mined by (6.1). Therefore, 'Harm is 1 -dimensional. 

(<=) If Au{n) = c,,-i n (u(n) - u(n - 1)) - c„,„ + i(m(« + 1) - u(n)) = for every n, then 

e„_i,„(u(n) - m(« - 1)) = c„,„ + i(m(« + 1) - w(n)) = a, 
for some fixed a (the amperage of a sourceless current). Then 



£(«) = Y c„_i >n (M(n) - u(n - l)) 2 = a 2 Y — !— < oo, (6.2) 

since u e 'Harm c *H S . Note that (6.2) implies m is bounded: £(w) = a 2«ez( M (") _M (' 1_ 1)) 
and £„>i(«(rt) ~~ M ( w ~~ 1)) = lim„^co u(ri) - u(Q). The function u is monotonic because it is 
harmonic, so the sum is absolutely convergent. □ 

Lemma 6.5. On (Z, c"), the energy kernel is given by 



v„(Jk) 



0, k < 0, 

1 < Jk < n, n > 0, 

1— r ' ' ' 



l_r»+i 



n \A) ( n \ — /7**M n ' — _ 

2(1-/-)' 



1-r 

and similarly for n < 0. Furthermore, the function w (n) — ar 1 " 1 , a := 7(l 'l r) , defines a 
monopole, and h(n) — sgn(«)(l — w c («)) defines an element of 'Harm. 

Proof. It is easy to check that Aw o (0) = 2c(a-ar) = 1, and that Aw () (n) = c"(ar" -ar"~ l ) + 
c" +l (ar" - ar" +l ) = for n + 0. The reader may check that &{w ) = 2 (\- r ) so tnat w o G He- 
The computations for v t and h are essentially the same. □ 

Remark 6.6. With h(n) = sgn(«)(l - w (n)) defined as in Lemma 6.5, the boundary term 
is 1. To compute this, use the exhaustion Gk = \—k, k], 

Y h(x)^(x) = 2c A '(ar t - a/" 1 ) = 2a(r- 1)- = 1, 

AGbd Ci r 

so that £ bdG /j£ = lim^oo 2bdG t /? £ = L 

Remark 6.7. As in Lemma 6.5, it is straightforward to check that w (n) = ar n \ a := tjz^, 
defines a monopole on the geometric half-integer model (Z+, c"). However, it is also easy 
to check by induction that 'Harm = for this model. 

Example 6.8 (Decomposition in D). In Remark 4. 14, we discussed the Hilbert space D and 
its inner product (u, v) := u{o)v{o) + (u, v)g. Since (Z+, c") and (Z, c") are both transient 
for c > 1 (but only the latter contains harmonic functions), it is interesting to consider Pd 1 
for these models (see Remark 4.13). The projections v = Pd 1 and u = 1 - v = ^d ^ on 
(Z, c") are given by 
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v(x) = 2-2a + ar' x < and u(x) = 2a - 1 - ar w , (6.3) 

where with a = =4r- and r = c _1 , and one can check u e At"; see Definition 4.11 and 
Lemma 4.12. In particular, Av = (1 - v )6o and Am = -u„5q (as usual, o = 0). Now 
consider the representative of w e M a given by 

w(x) = (2- <!)*[_„, 0] + (l + 2a(r |A| - c))* [li0o) . (6.4) 

A straightforward computation shows that w — v + h with h e HD . 

The function v = Pd 1 was computed for (Z, c") in (6.3) by using the formula £(w) = 
u a — Ug, from Lemma 4.15, where u := ^d 1 = 1 - v and u a = u{o). For a general 
network (G,c), this formula implies that (u ,&(u)) lies on a parabola with u e [0, 1) and 
maximum at (i, |). From (6.3), it is clear that the network (Z, c") provides an example of 
how u„ = 1 - T^-j- can take any value in [0, 1). Note that c = 1 corresponds to fi(«) = 0, 
which is the recurrent case. 

Example 6.9 (Star networks). Let (S,„,c") be a network constructed by conjoining m 
copies of (Z + , c") by identifying the origins of each; let o be the common origin. 

In [JP09a], we explore the boundary bdG in more detail. The idea is that the boundary 
term is nontrivial precisely when bdG + 0. The presence of a monopole indicates that 
bdG contains at least one point; see Theorem 4.6. If "Harm + 0, then there are at least two 
boundary points; see Lemma 5.5 and Corollary 4.4. 

Example 6.9 shows how to construct a network which has a boundary with cardinality 
m. Note that these boundary points can be distinguished by monopoles as in (6.4). This 
monopole acts as a sort of indicator function for the corresponding boundary point +oo; 
when we pass to TYg, w is supported only on the positive half of Z. We return to the 
general case (G, c) in [JP09a]. 

Example 6.10 (Independence of exhaustion). The following elementary example shows 
how the boundary term may not be independent of exhaustion if one does not restrict to 
u € M; this example is adapted from [KY89, Ex. 3.12]. Take the nonnegative integers 
{0, 1, 2, 3, ... } with unit conductance edges connecting nearest neighbours, as depicted. 

1111 
• • • • ■ ■ ■ 

12 3 

Define a function u on the vertices by 

f l/k, if n — 2 k , 
m(0) = and u(ri) = u(n — 1) + < 

ll/n, otherwise, 

so that the increment between neighbours is either - or \ — r^— . The energy is estimated 

& n k log 2 n & - 7 

OO ^ CO j 

which proves that u e 7Yg. Note, however, that u is unbounded, so that u t At. Let {Gjtj^lj 
be an exhaustion by the sets Gk = [0,2*], so that bdG^ = (2*). Then the boundary sum 
contains only one term and we have 
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\ 7=1 7 J 



bdG 



--log 2. 



1 On the other hand, let \G' k }™ =l be an exhaustion by the sets G' k = [0, 3 ], so that bdGjJ, 

2 {3 A }. Then the boundary sum still contains only one term but now we have 



Z 

bdG 



u(3 k )^(3 k ) < 



I V 

Z^ 



W=i 



1 

3 1 



-> 0. 
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